Detecting phonons and persistent currents in toroidal 
Bose-Einstein condensates by means of pattern formation 
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We theoretically investigate the dynamic properties of a Bose-Einstein condensate in a toroidal 
trap. A periodic modulation of the transverse confinement is shown to produce a density pattern 
due to parametric amplification of phonon pairs. By imaging the density distribution after free 
expansion one obtains i) a precise determination of the Bogoliubov spectrum and ii) a sensitive 
detection of quantized circulation in the torus. The parametric amplification is also sensitive to 
thermal and quantum fluctuations. 
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Bosc-Einstcin condensates have recently been obtained 
with ultracold gases in a ring-shaped geometry [1] . Other 
groups are proposing different techniques to get toroidal 
condensates [J. The aim is to create a system in which 
fundamental properties, like quantized circulation and 
persistent currents, matter-wave interference, propaga- 
tion of sound waves and solitons, can be observed in a 
clean and controllable way. An important issue concerns 
the feasibility of high-sensitivity rotation sensors. 

In this work we show that key properties of conden- 
sates in toroidal traps can be measured by means of para- 
metric amplification. This corresponds to an exponential 
growth of some excited modes of the system induced by 
a periodic modulation of an external parameter We 
consider the modulation of the transverse confinement, 
which is shown to drive a spatially periodic pattern in 
both density and velocity distributions as a consequence 
of the amplification of pairs of counter-rotating Bogoli- 
ubov phonons. If the trap is switched off, this pattern 
produces a peculiar flower-like density distribution of the 
freely expanding gas. The number of "petals" and their 
shape provide a sensitive measure of the excitation spec- 
trum and the superfluid rotation of the condensate. 

We perform numerical simulations by integrating the 
Gross-Pitaevskii equation [J] for N bosonic atoms of mass 
M, confined in an external potential V ex t- 



ihdt'ip = 



2M 



V 2 + V a 



(i) 



The mean-field coupling constant is given by g — 
AiiTi 2 a/M, where the s-wave scattering length a is as- 
sumed to be positive. The order parameter of the con- 
densate, ip(r,t), is normalized to J dr\ifj\ 2 — N and may 
by written as i)j — n 1 / 2 exp(i5), where n is the conden- 
sate density and the phase S is related to the superfluid 
velocity by v = (h/M)VS. We consider a condensate of 
N = 2 x 10 5 atoms of 87 Rb confined in a torus of length 



L = 2nR = 100 fim by an axially symmetric potential 
which has a minimum, V ex t — 0, at z — and r± = R. 
The trap is harmonic and isotropic in the (z,r^)-plane 
around this minimum, with frequency u± = 2tt x 1 kHz. 
With this choice the transverse width of the condensate, 
ro, is significantly smaller than the radius of the torus. 
This implies that curvature effects are almost negligible 
in the ground state and in the in-trap dynamics and one 
can replace the torus of radius R with a cylinder of length 
L with periodic boundary conditions [5]. Curvature ef- 
fects are instead important during the free expansion of 
the condensate, and therefore the full toroidal geome- 
try is used for simulating the dynamics after the release 
from the trap. The choice made for the geometry and the 
parameters is intended to simulate feasible experiments, 
but the effects we are going to show can be obtained with 
different types of traps and in a wide range of parameters. 

The process that we want to study has three steps: i) 
the condensate is initially prepared in the torus; ii) the 
transverse harmonic potential is periodically modulated 
for a time i mo d; iii) both the trap and the modulation are 
switched off and the gas freely expands for a time t exp . 

We first solve the stationary GP equation and the Bo- 
goliubov equations obtained by linearizing Eq. |T]) [J] for 
a cylindrical condensate, by using the same numerical 
techniques as in We obtain the stationary solution 
^(r) and its excitation spectrum, namely the eigenfre- 
quencies u>i and quasiparticle amplitudes Uj(r) and Vi(r). 
The condensate at step (i) is assumed to be at equi- 
librium. Quantum and/or thermal fluctuations are in- 
cluded by means of the Wigner representation of bosonic 
fields In practice, the function ip, input of step 

(ii), is taken of the form ip — ipo + J2ii c i u i + c i v t\i 
where the sum extends over a wide set of Bogoliubov 
states including those which are expected to be rele- 
vant for the subsequent parametric amplification. The 
Wigner representation enters through the choice of the 
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to occur via the parametric instability of Bogoliubov 
phonons, with wavevector k, propagating along the torus 
in opposite directions. The mechanism is basically the 
same already described in 11, Ti| for elongated conden- 



FIG. 1: a) Density of a condensate in the toroidal trap af- 
ter a transverse modulation with f2 = 0.6 ui±, A = 0.1, and 
tmod = 130 ms, corresponding to a maximum of visibility of 
the density pattern; the density exhibits 8 oscillations and the 
relative density variation is about 20%. b) Tangential current 
for the same condensate, c) Density after 7 ms of free expan- 
sion, d-e-f) Same as above but for t mod = 130.8 ms, when 
the pattern visibility is maximum in the current distribution. 
Each box is 45 /im per side. Densities and currents are in- 
tegrated along the axis perpendicular to the torus. White 
(black) means high (zero) density, while in frames (b) and (e) 
white and black correspond to currents in opposite directions. 



complex coefficients a, which are taken to have a ran- 
domly distributed phase, with constant probability in 
[0, 2ir), and modulus, with Gaussian distribution of width 
1/2 + [1 - expinui/ksT)]- 1 . 

To simulate the in-trap dynamics at step (ii) we nu- 
merically integrate Eq. ([1]) [8J. The external modula- 
tion is included by multiplying V cx t by the factor f(t) = 
[1 + s(t)A cos(flt)] , where A is the amplitude of the mod- 
ulation and s(t) is an appropriate switching function. A 
collective motion which could be easily excited by this 
perturbation is the radial breathing mode of frequency 
2w_l @- However, the processes that are of interest here 
are those involving longitudinal phonons at frequency 
lower than 2uj± . We also want that the transverse width 
and the density of the condensate oscillate in-phase with 
the external modulation. To this aim, the possible cou- 
pling with the breathing mode has to be suppressed. This 
can be obtained by choosing s(t) to be sufficiently adia- 
batic. In the following we will use an amplitude A = 0.1, 
for which we find that a switching time of about 20 ms 
is adequate in a wide range of Q. The entire modulation 
time, t m od, is instead chosen to be longer than 100 ms. 

The transverse modulation is found to induce a longi- 
tudinal pattern formation [10| corresponding to a stand- 
ing wave along the torus. An example is shown in 
Fig. QJa) where the in-trap density distribution is plotted 
for a modulation of frequency Q — 0.6a->^ and duration 
tmod = 130 ms. In this case the density oscillations along 
the torus have 8 maxima and correspond to a relative 
density variation of about 20%. This pattern is found 



sates in a modulated ID optical lattice and it is analogous 
to the spontaneous pattern formation discussed in [l3| in 
different geometries. It has also interesting similarities 
with the phenomenon of Faraday's instability 13, 14 | for 
classical fluids in annular resonators [lj| . 

The spatial periodicity of the standing wave in the 
torus is fixed by the wavelength of the most rapidly grow- 
ing modes, determined by the Bogoliubov dispersion re- 
lation Lo(k) and the boundary conditions. The paramet- 
ric amplification spontaneously selects modes with fre- 
quency close to the resonance condition uj(k) = fi/2 and 
wave vector k = 2-7rm/L, where the integer m represents 
the azimuthal angular momentum of the excitation and 
2m is the number of nodes in the density pattern. The 
position of the nodes changes randomly at each realiza- 
tion, being related to the phase of the initial fluctuations 
which are parametrically amplified. 

At resonance the ±fc components of the Fourier trans- 
form of the order parameter along the torus, ip(k,r±), 
grow exponentially during the modulation. Since the 
GP equation JT]) contains a nonlinear mean-field term, 
the growth is limited by the energy transfer to non- 
resonant modes due to mode-mixing processes. Thus 
the ±fc components eventually saturates around a max- 
imum value for long times. An example is given in 
the upper panel of Fig. ® where we show the quan- 
tity Pfc = 2ir J dr± r±\ip(k, r±)\ 2 as a function of t m od 
for Q = 0.6lu±. In the lower panel we plot the maxi- 
mum value achieved by in case of a modulation with 
tmod up to 200 ms, as a function of ft. The position of 
each peak coincides with twice the frequency of the cor- 
responding m-phonon, as calculated by solving the Bo- 
goliubov equations. The figure suggests the possibility to 
use the parametric amplification for spectroscopy: Each 
time a density pattern is observed, the wavevector k is 
simply obtained by counting the number of oscillations 
in the torus, while the frequency u(k) is just fl/2 [l6j |. 

The standing density wave is phase-locked with the ex- 
ternal modulation: the amplitude of density oscillations 
is maximum when the width of the condensate is cither 
maximum or minimum, and vanishes when the width is 
the same of the unperturbed condensate. This also causes 
the fast oscillation of the quantity Pfc(t) in Fig. [2] The 
plots in Fig. QJa) and (d) correspond to a maximum and 
a minimum of Pfe(t), respectively; the density fluctua- 
tions in Fig. [ljd) are random and small. The standing 
density wave is accompanied by a pattern in the velocity 
distribution. This is easily seen by plotting the tangen- 
tial current J = (ih/2M)[ip(->p*y - (ip)'ip*], as we did in 
Fig. [TJb) and (e) : the velocity distribution exhibits small 
random fluctuations in (b) and a clean periodic pattern 
in (e). The density and velocity distributions oscillate 
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FIG. 2: Top: Amplitude of the ±fc components of the Fourier 
transform of the order parameter along the torus, where the 
wavevector k satisfies the resonance condition co(k) — $7/2. 
The amplitude is plotted as a function of t mo d for $7 = 0.6u>± 
and A = 0.1. Bottom: maximum value of in the interval 
< tmod < 200 ms as a function of $7 (in units of uj±). 
Solid squares on the horizontal axis correspond to twice the 
eigenfrequencies u)(k) of the Bogoliubov equations, with k = 
2nm/L and m integer. 



exactly out-of-phase in time. 

The occurrence of a pattern in the velocity field has 
spectacular consequences in the expansion after the re- 
lease from the trap. In fact, the presence of ±fc-phonons 
gives rise to interference fringes of atoms expanding in 
preferred directions, similar to those observed in the ex- 
pansion of an elongated condensate with Bragg excited 
phonons [T7[ ■ In toroidal geometry these fringes produce 
a flower-like structure with m "petals" reflecting the pe- 
riodicity of the initial velocities. An example is given in 



Fig. [ljc) and (f) where we show what happens to the 
condensates of Fig. QJa) and (d), respectively, when the 
trap is switched off and the gas is imaged after 7 ms of 
free expansion. The figure shows that, by starting the 
expansion when the pattern has maximum amplitude in 
velocity field, one obtains very clean interference fringes. 
Since these "petals" are much more visible than the in- 
trap density oscillations, the expansion significantly en- 
hances the sensitivity of the spectroscopic measurement. 

The dynamics in Fig.[ljc) and (f) is simulated by solv- 
ing Eq. U) without the external potential and neglect- 
ing the mean-field interaction. From the experimental 
viewpoint, this situation corresponds to either tuning the 
scattering length a to zero by means of a Feshbach reso- 
nance just before the expansion or starting from a torus 
strongly squeezed in the z-direction. In the latter case, 
due to the fast expansion along z, the density decreases 
much faster than the typical timescale for the transverse 
expansion, thus suppressing mean-field effects in the ra- 
dial motion. However the choice to neglect the mean- 
field interaction is not necessary for obtaining a distinct 



FIG. 3: Density of the expanded condensate as in Fig. QJf). 
The expansion starts again from the condensate in Fig. [TJd) , 
but here the mean-field interaction is included in the simula- 
tion of the expansion, while in Fig. [TJf ) it is neglected. The 
expansion time is t Exp — 1.8 ms. The box is 45 /im per side. 



flower-like structure. This can be seen in Fig. [3] where we 
plot the density after free expansion including the mean- 
field interaction for a condensate prepared as in Fig.[TJd); 
the petals are still visible also in this case. The main dif- 
ference with Fig. QJf) is the timescale of the expansion, 
which is faster in the presence of the repulsive mean-field 
interaction. A rough estimate of the ratio of the two 
timescales, which holds at short times, is easily obtained 
by comparing the expansion of a cylindrical condensate 
in these two cases: i) the scaling behavior of the full GP 
equation 0] and ii) the dispersion of a noninteracting 
wavepacket with the same initial shape. Using the pa- 
rameters of our condensate, the interacting gas turns out 
to expand faster by a factor 7 ~ 3.3. In a toroidal geome- 
try this factor is expected to be slightly larger. The ratio 
between the two expansion times in Fig. [T] and Fig. [3] 
has been taken to be ~ 4 and the size of the expanding 
condensate and of the petals is indeed similar. 

The process of pattern formation is significantly af- 
fected by the presence of quantized circulation. If the 
condensate is initially rotating with angular momentum 
L z = nh per particle, where k is the quantum of circu- 
lation, then the external modulation produces a pattern 
which rotates along the torus at the same angular veloc- 
ity of the condensate as a result of the frequency shift of 
counter-rotating phonons [l8| . More strikingly, the inter- 
ference structure observed in the free expansion exhibits a 
significant misalignment of opposite petals, proportional 
to k. This is evident in Fig. |4] where we show the same 
expanded condensate but for different values of k. The 
axis of each petal can be approximated with a rectilinear 
segment tangent to a circle of radius i?o, so that 2R a is 
the distance between the axis of opposite petals. The 
radius Rq is of the order of the size of a k- vortex core in 
the expanding condensate. Its time evolution is well re- 
produced by the law Ro(t) = Rna + /3nt/R. The second 
term, with of order 1, is consistent with the asymptotic 
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FIG. 4: Interference patterns in the expansion of a paramet- 
rically excited condensate in the presence of quantized circu- 
lation (k = 1,5, 10 from (a) to (c)). The superfiuid rotation 
causes a misalignment of the fringes proportional to k. In all 
plots t exp = 7 ms and t mo d = 180.8 ms. Each box is 55 /im 
per side. 



expression for the vortex profile that can be derived ana- 
lytically in the case of free expansion as discussed in [l9| . 
The quantity na in the first term is of the order of the 
ratio between the rotational velocity of the condensate 
in the torus and the velocity at which the condensate ex- 
pands radially; this ratio is not ~ (hK/MR)/(u>±ro/^f)- 
Fig. |4] shows that the parametric amplification process 
can be used as a sensitive probe of rotations in the torus 
[20| , which works down to a few quanta of circulation and 
is complementary to the techniques discussed in [l9j |. 

All results shown above have been obtained assum- 
ing the condensate to be at T = 0. In this case, the 
fluctuations which are parametrically amplified are those 
due to the small quantum depletion of the condensate. 
However the dynamics of the amplification is greatly in- 
dependent of the type and origin of the initial fluctu- 
ations. What really matters is the amount of effective 
seed (i.e., a specific quadrature of counter-propagating 
Bogoliubov excitations [12|]), which determines the time 
tmod needed to obtain a visible density pattern. In an ex- 
periment the imperfections in loading the condensate in 
the trap and thermal fluctuations can also contribute to 
the seed. In particular, one expects that thermal fluctu- 
ations become dominant for temperature T 3> Ml/kB- 
We have performed simulations in this regime finding 
that the time i mo d needed to achieve the maximum pat- 
tern visibility, defined as in Fig. [21 scales approximately 
as (l/n)ln(Q/T). This is consistent with the fact that 
the growth rate and the amount of seed are expected 
to be proportional to f2 and T/Q, respectively 12j. For 
the parameters of Fig. [2] and for T » SO/fc^ a factor 2 
in temperature implies a change of about 10 ms in the 
timescale of the amplification. This suggests that the 
visibility of the pattern as a function of $1 could provide 
a way to measure T in a regime of temperatures where 
standard thermometric methods are not applicable. 

In conclusion, we have shown that a periodic mod- 
ulation of the confining potential of a toroidal conden- 



sate induces a spontaneous pattern formation through 
the parametric amplification of counter-rotating Bogoli- 
ubov excitations. This can be viewed as a quantum ver- 
sion of Faraday's instability for classical fluids in annular 
resonators. The occurrence of this pattern in both den- 
sity and velocity distributions provides a tool for mea- 
suring fundamental properties of the condensate, such as 
the excitation spectrum, the amount of thermal and/or 
quantum fluctuations and the presence of quantized cir- 
culation and persistent currents. 
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